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We study full counting statistics of coherent electron transport through multi-terminal interacting 
quantum-dots under a finite magnetic field. Microscopic reversibility leads to the symmetry of the 
cumulant generating function, which generalizes the fluctuation theorem in the context of quantum 
transport. Using this symmetry, we derive the Onsager-Casimir relation in the linear transport 
regime and universal relations among nonlinear transport coefficients. 
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Full counting statistics (FCS) has become an active 
topic in the mesoscopic physics 0, HI EL H, IS S B 



ULfl, 11 



12j . FCS addresses the probability distribu- 
tions of charges transmitted during a measurement time, 
t. It can characterize the statistical properties of quan- 
tum transport in the far-from-equilibrium regime. Since 
the first paper by Levitov and Lesovik many the- 
ories have clarified various aspects of the distributions 
0, [H, H, 0, H, 0| ■ Recently, experiments have been con- 
ducted to measure third current cumulants [l(| and the 
distributions [H|,[l2j]. However, FCS has never been ap- 
plied to exploration of general aspects, such as nonequi- 
librium thermodynamic structures in coherent electron 
transport. In this paper, we discuss these general aspects 
by studying symmetries in FCS, which is valid beyond a 
linear response regime. 

Our argument is based on the microscopic reversibil- 
ity, and is related to the steady state fluctuation the- 
orem (FT) in nonequilibrium statistical mechanics [HI, 



141 . Il5l . |lq |. FT is an important theory that holds 
even in the far-from-equilibrium regime. It is written 
as ]ixa T ^ 00 lix[P(AS)/P(-AS)]/T = I E , where P{AS) is 
the probability of entropy AS = Iet, produced during 
time r. This expression quantifies the probability of neg- 
ative entropy, which can be finite for a short interval of 
time in small systems, as demonstrated in the context of 
the second law violation in colloidal particle experiments 
Remarkably, FT can reproduce Onsager's recipro- 
cal relations and the Kubo formula lH, 17 1 and predicts 
properties in the far-from-equilibrium regime [17j |. 

Recently FT was studied with regard to classical meso- 
scopic electron transport, i.e., in the sequential tunneling 
regime in th e q uantum-dots, with Markovian approxi- 
mations [l?], [3- These works highlighted the relation 
between FT and FCS. In this paper, we study the gen- 
eral relation between FT and FCS with respect to coher- 
ent electron transports in generic situations. We consider 
a multi-terminal interacting electron cavity under a fi- 
nite magnetic field. The validity of FT has not yet been 
established in this regime. In electron transport, two 



thermodynamic forces -the thermal gradient and the bias 
voltage- produce entropy. Thus, we introduce a cumulant 
generating function to obtain heat and charge transfer, 
and exactly derive its new symmetry. The symmetry 
leads to a quantum version of FT and measurable uni- 
versal relations among nonlinear transport coefficients. 
They are extensions of the Onsager-Casimir relation [l9j |. 

Model. — We consider a mesoscopic cavity connected 
to m electron reservoirs. The total Hamiltonian consists 
of the reservoirs H r (r = 1, • • ■ , m), cavity Hd, interaction 
i?int, and tunneling Ht 



H — E™ 1 Hr + H c i 



H int + H T . (1) 

= Etfff 4a d 3<y Where 



The cavity is described as Hd 

die annihilates an electron with spin a at site i. We 
assume the Coulomb interaction in the cavity H xu ^ = 
Hijaa' Uiaja' d\ a d i°' d ^ A Tnc Hamiltonian for the 
reserver r is H r = J2ka e rka>l kcr a r ka, where a r ha annihi- 
lates an electron with spin a and the wave vector k. The 
tunneling between the reservoirs and the cavity is de- 
scribed as H T = J2 r kia l rki d\ a a rk a + H.c.. When a mag- 
netic field B is applied, the elements of the hopping and 
tunneling matrix acquire the phases ty = \tij\exp(i(f>ij) 
and t r ki = \t r ki \ exp(i <pri). The phases are odd functions 
of the magnetic field <j>(—B) = —(f)(B). 

Our calculations follow the standard procedure with a 
perturbation series for H lvA . Throughout this paper, we 
use h = fee = c = 1. The density matrix at the initial 
time t = — t/2 is assumed to be of the product form po = 
Y\ s p s , (s = 1, • • • , m, d), where p s is the equilibrium dis- 
tribution at temperature T s — \/f3 s and chemical poten- 
tial p, s ; p s =exp[- f3 s (H s -p s N S )]/Tr exp[- (3 S (H s -p, s N s )]. 
The operator N s is the number operator in reservoir 
N r = J2ka a lka a rka, or the cavity N d = E w i4. 
The charge and heat current operators are defined as 
I cs = i [N s , Ht] and Its = i [H s , Ht]- The expression for 
currents I cr and I^r are given by 



Ihr = — J2 
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H.c. , 
+ H.c. 
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Cumulant generating function. — We introduce the 
characteristic function (CF) for the transmitted charge 
and heat during time r, q as — Jll^dt la s(t) (a = c,h): 

Z({ Xas };B)= (V^e iH -V 2 e- iHT V^), V = \[,V» (2) 

where V s =exp[— i(xhsHs+XcsN s )/2], This contains the 
counting fields for charge and heat current Xcs and Xhs, 
respectively. The symbol (...) denotes an average over the 
initial state. Equation ([2]) can be rewritten in the familiar 
form, where the counting fields play roles of fictitious 
gauge fields [H, 0, since Vj H T V r makes the gauge 
fields in the elements of the tunneling matrix. Note that 
Eq. © satisfies the normalization condition Z({0}; B) = 
1. We use the CF and introduce the cumulant generating 
function (CGF) defined at the stationary state as 



where g(cu) is the Green function which has site, spin, 
and Keldysh indices. The normalization condition 
J-q({0}; B) = must be satisfied. The matrix elements at 
(io~,jo~') for the inverse of the Green function are given 
in the Keldysh space as 



[9 'L^^^ 



Ujf 3 ~Y, r i ^rij{u)f 3 - f/ d , 



for a = a' and for other cases. The self-energy and the 
Pauli matrix f 3 are 



H r ij(u>) 



1 

-1 



In the wide-band limit T r ij = 2w^2 k t r ki trkj^( L0 ~ £ rk) |22l |. 
the causal and anti-causal components are S • • (w) 
I [1/2 — /+(w)] with the Fermi/hole distribution 
•^({Xar}; B) = linv^oo lnZ({x Q r}; B)/t. (3) function /±(a/j = l/{exp[±/3 r (w - fi r )\ + 1}. The lesser 

and greater components depend on reservoir's counting 
fields as £±J( u) = ±ir rij f^(uj)exp{±i(xhrU + Xcr)}- 
The wide-band limit is not critical to our argument as 
long as the reservoirs form continuum energy spectra, 
which cover the entire energy range of the cavity. The 
term f/d has an infinitesimal contribution, depending on 
the counting fields and the initial state of the cavity. This 
term is crucial in the causality 2l|, but negligible com- 
pared to the self-energy terms of the reservoirs. This is 
plausible because the steady state should not depend on 
the initial state of the cavity. In addition, we consider 
the rotation of the Green function in the Keldysh space 
ivgicrj<T'(w)R, where the operator R is defined as 



The CGF generates cumulants from the derivatives with 
respect to the counting fields. The first and second 
derivatives generate the average current between the ter- 
minal and the cavity, and a symmetrized current corre- 
lation expressed as 

«/ cl » = ™ = lim 



dix 



cl 



T — >oc T 



(/J cl J c2 )\ = d 2 T({0};B) = Um ({qcl,qc2})-2(q cl )(q c2 ) 



diXddiXc 
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Microscopic reversibility. — So far several symmetries 
in CGF are known [9(. For instance, CGF is a 27r-periodic 
function of Xcs, which is a consequence of the discrete- 
ness of the charge. We take the time reversal symmetry in 
CGF into consideration. Let be the time reversal oper- 
ator, which evaluates the complex conjugate for complex 
numbers and reverses the spin operator [20j ]. This time 
reversal operator satisfies the equations; 0z0 _1 = — i 
and (n\0\n') = (h'\e0^e- 1 \h) $2%. Here O is an oper- 
ator and |n) = 0|n). Calculations with these equations 
and use of the definition of operator V lead to the relation 



R = exp[i(xhm.V + Xcm)T 3 /2] 



(6) 



This transformation does not change Eq. ([5]), instead it 
only replaces Xar by Xar — Xam in the Green function. 
Therefore, we conclude that J-q depends only on the dif- 
ferences between the reservoirs' counting fields. 

We now consider the interaction part T- ln t based on the 
linked cluster expansions. It is formally written as 



Z({ Xas }; B) = Z({- Xas + iA as }; -B) 



(4) 



Jint = lim In 



-isr 



JSj 



where A cs =f3 s p s and Ahs — ~Ps- This equality is critical 
in obtaining our central result ((9]). 

When stationary, no extra charge and heat accumulate 
inside the cavity. This implies that CGF depends only on 
the differences between the reservoirs' counting fields and 
is independent of the cavity's counting field and initial 
state. Below, we present its proof using the Schwinger- 
Keldysh approach [IBED]. 

The entire CGF expression consists of a noninteracting 
part J-q and an interacting part jT int . First we consider 
the noninteracting part. After a straightforward calcula- 
tion, the analytical expression of J-q is obtained as 



Here the functions Sj and S^L are given as 



It (7) 
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dtdt'Jivitf^givjvi^t'^Jjvit'), 



-r/2 



r/2 



dtH iv 



r/2 



SJ* a± (t)' 8J ia ±{t) 



1 f°° 

FodXaa}] B) = -— / dui In detg(w) + const. 

27TJ-OC 



(5) 



where the Grassmann source field in the Keldysh space 
J%a = t (Jicr+, Jia~) is used. The function i?int [•--] is ob- 
tained by substituting the derivatives of the Grassmann 
numbers for the fermion operators in the Hamiltonian 
-Hint- We consider these functions in the Fourier space, 
and transform the Green function and Grassmann fields 
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using the rotation operator R. Then we readily find 
that ^int is expressed by the simple replacement of the 
Green function in Sj by the inverse Fourier transform of 
R^9iajty' (u)R- In this mathematical structure, the charge 
and energy conservation in the interaction process is cru- 
cial. To check it, let us consider the calculation for one 
diagram that contributes to the second order in J-- mt 



- \ T 3 T 3 

4 Z^i s s s ' 



'duj dv dv 1 



X 9i'a,ia 



2n 
A v ')9 



^(i/+w),(8) 



s s 



where the summations are performed over the site, spin, 
and Keldysh indices s, s' = ±. Eq. §8§ does not change 
on the rotation Wgi a j a > (ui)R, since diagrammatically the 
energy and charge are conserved. Therefore we find that 
Tint depends only on the difference of the reservoirs' 
counting fields. From Eq. ((31) , we obtain the symmetry 
in the interacting electron transport as 



H{x ar }; B) = T({- Xar + tA ar }; -B), 



(9) 



where Xar = Xc 



An 
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Xam and A ar represents the affinity 
We can check that CGFs in interacting 



quantum-dots in Refs.@, 0, i, @] satisfy Eq. © 

Fluctuation theorem. — Eq. © can be regarded as 
a quantum fluctuation theorem (FT). Let us consider a 
two-terminal setup and define the entropy produced as 



AS = Ac\q c \+Ahiqhi- 



(10) 



We fix the counting fields as Xc\A~^ — XhiA^l = X- 
The entropy production is obtained from the derivative 
of CGF with respective to x- Asymptotic form of the 
probability distribution can be obtained from the inverse 
Fourier transform of the CF with the variable x- Saddle- 
point analysis with the symmetry ([9]) yields the relation 



lim In 



P(AS; B) 



P(-AS; -B) 



It = I E 



(11) 



This formula generalizes FT in the quantum regime un- 
der a finite magnetic field. At a uniform temperature, 
the entropy production is proportional to the charge cur- 
rent. In this case, Eq. (fTTj) quantifies the probability of 
back flow of charge currents. 

General relations in nonlinear transport regime. — 
We show that the symmetry © predicts general rela- 
tions among nonlinear transport coefficients. We con- 
sider the situation that only the chemical potentials 
/ii and /i2 are varied at a uniform temperature /3 , 
and the charge currents are measured at terminals 1 
and 2. We then compute cumulants of the currents, 
((Iciloi)) = d k ^T({0};B)/d(zXci) k ^^Xc2) k2 . The 
nonlinear transport coefficient is defined in the expan- 
sion of the cumulants with respect to the affinities as 



L 



(B) 



d< 



! «4i4 2 » 



dAHdAl 
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FIG. 1: (Color Online) The flux dependence of nonlinear 
transport coefficients. Parameters: £ = /i, t = 10F, and /3 = F. 
Eas. (ri5jl and (Tl7|) are satisfied. is always zero. 



The affinities A C j(j = 1, 2) are written as A C j = P(pj — fjf) 
with a definite chemical potential \x for terminals 3 to m. 
We symmetrize the coefficients and CGF as 



kik-2 



= L 



(B) ± L%%(-B) 



-}) = T({XcrhB)±T({ X crY,-B) 



From Eq.©, we obtain the equality fF±{{Xcr}) = 
±fF±{{ — Xcr + iArr})- Note that CGF is a function of 
A cr as well as Xcr [171 ]. By taking the derivatives with 
respect to the affinities and counting fields for both sides 
of the equality of J-±, we obtain the general relations 
among the coefficients as 



n 2 =0 



(-1) 



n n i^\* 2 +£ 2 i± ,(i3) 



J Cl-Ill f. 2 



where n = ni+n2+ki+k2- With the equalities ± = 0, 
the relations (|13j) can be further simplified. We consider 
equations (fl3|) which satisfy Af — ki+k2+£ Equations 
for Af — 2 reproduce the linear response results such as 
Kubo formula and Onsager-Casimir relations [23} 



L%(-B), L\l{B)=L\l{-B), (14) 



and Lq\(B) = Lq\(—B). Relations beyond linear re- 
sponse regime can be obtained for Af > 3. We list some 
of the relations for Af = 3 as 
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(15) 
(16) 
(17) 
(18) 



In mesoscopic experiments, large bias voltages can easily 
produce finite higher order coefficients, and the Onsager- 
Casimir relations can be violated 25, 2(| 27|]. Eqs. (fT5|) - 
(fTHI) demonstrate that beyond the Onsager relation, uni- 
versal relations exist in the nonlinear transport regime. 
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These nontrivial relations rely solely on the microscopic 
reversibility and thus are insensitive to the setup details. 

Three-terminal Aharonov-Bohm interferometer. — The 
simplest setup to demonstrate these relations would be a 
three-terminal Aharonov-Bohm (AB) ring with a three- 
fold symmetry. We consider a ring consisting of three 
nonintcracting quantum-dots, each of which connects to 
a reservoir, as shown in the inset in Fig. [T] The Hamil- 
tonian of the ring is given as Hd = J2 a Si=i £ d\ a d icr — 
te^^dh^dicr +H.c. {d& a = d\a). The tunnel coupling is 
T r ij = T5ijd r i. The explicit form of Eq. ([5]) is given by 



J(Xj-Xk) . 



1) 



1 f°° ( 
T = - do; In 1+]T //(w)/ fc ->)(e 

J -°° L j,k=i 

x[% vcn (uj) - y^T^ddM ejiujl - 2 //(w))] L 
£=1 J 

where £jfe£ is the totally antisymmetric tensor. The 



transmission probability T e% 



a/odd 



are written as T even = 

- f3, 



P(l/4 + P + z 2 - 2izcoscj))/A and T odd = Psm<p/A, 
where z=(u-e)/T, t = t/T and A = (z 2 +l/4) 3 +P(3/8- 
6z 4 )+t 3 z[(4z 4 -3)-12t 2 ]cos</>+9t 4 (z 2 +l/4)+4f 6 cos 2 ./). 
Figure [1] shows the linear conductance multiplied by 

[J + , as well as the nonlinear coefficients in Eqs. (fl"5]) and 
(jTTJ) . It is clear that the relations (fTSf and (fT"7|) are sat- 
isfied. The overall structures depend on the temperature 
regions. 

Two-terminal case. — Two-terminal geometry is a 
common setup used in experiments. Eqs. (fl~5"|) - (fl"8l) 
are still satisfied, but the double script notations for 
the transport coefficients are replaced by single script; 



L 



fel k 2 



(_i)(fa+«i) L fci 



fcl+fe2 



In general, all cumulants 



of charge currents are symmetric in the magnetic field in 



nonmteractmg systems 
cannot be measured 



24. 25 



26J. Hence, Eqs.p7 | -(fT8 l) 

in noncentrosymmetric 
is generally finite 



However 

interacting systems, the coefficient L\ 
[28j |. In this case, Eq. (|18p predicts a finite third current 
cumulant even in the equilibrium. We confirmed a finite 
value of the coefficient L\ _ and Lq _ which is asymmet- 
ric in the magnetic field, in the model of a double dot AB 
interferometer with noncentrosymmetric geometry where 
one of the dots is capacitively coupled to an additional 
gate electrode. We expect that the relation (fi~8|) can be 
observed in the setup of the experiments [27| . 

Summary. — We derived the new symmetry in full 
counting statistics from the microscopic reversibility. 
This symmetry can be regarded as the quantum version 
of the fluctuation theorem under a magnetic filed. When 
there is no magnetic field, this reproduces the usual fluc- 
tuation theorem without quantum corrections. The sym- 
metry leads to relations among nonlinear transport coeffi- 
cients. Eqs.(fl~5 |) -(fT8 |) would be measurable in present-day 
experiments. We hope this work will motivate further 
experimental and theoretical work on universal relations 
among nonlinear transport coefficients. 
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